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Highly-deformable materials, from synthetic hydrogels to biological tissues, are becoming increas¬ 
ingly important from both fundamental and practical perspectives. Their mechanical behaviors, 
in particular the dynamics of crack propagation during failure, are not yet fully understood. Here 
we propose a theoretical framework for the dynamic fracture of highly-deformable materials, in 
which the effects of a dynamic crack are treated with respect to the nonlinearly deformed (pre¬ 
stressed/strained), non-cracked, state of the material. Within this framework, we derive analytic 
and semi-analytic solutions for the near-tip deformation fields and energy release rates of dynamic 
cracks propagating in incompressible neo-Hookean solids under biaxial and uniaxial loading. We 
show that moderately large pre-stressing has a marked effect on the stress fields surrounding a 
crack’s tip. We verify these predictions by performing extensive experiments on the fracture of soft 
brittle elastomers over a range of loading levels and propagation velocities, showing that the newly 
developed framework offers significantly better approximations to the measurements than standard 
approaches at moderately large levels of external loadings and high propagation velocities. This 
framework should be relevant to the failure analysis of soft and tough, yet brittle, materials. 


I. BACKGROUND AND MOTIVATION 

Material failure is mediated by the propagation of 
cracks, which are spatially-extended dissipative defects 
that concentrate large deformation and stresses near their 
tips. The latter is an essential physical property of cracks 
that highlights their basic role in material failure: cracks 
strongly amplify remotely applied stresses near their tips, 
precipitating fracture as they propagate. 

The classic approach to the fracture of brittle materi¬ 
als, Linear Elastic Fracture Mechanics (LEFM) [IHS], is 
a perturbative approach that treats all relevant physical 
quantities to linear order in the elastic deformation with 
respect to the undeformed state of the material. Here, all 
nonlinearities and dissipation are assumed to be confined 
to a microscopically small region near the crack tip, and 
are neglected. 

LEFM does not treat separately the effects of the ex¬ 
ternally applied loading (pre-stress/strain) and the pres¬ 
ence of the crack itself, due to the linearity assumption. 
This physical picture makes sense as long the materials 
of interest fail when the applied forces are small and the 
deformation is indeed linear elastic (i.e. infinitesimal) ev¬ 
erywhere except for the immediate vicinity of the crack 
tip. This has always been the case in “traditional” brittle 
materials such as glasses, ceramics and brittle polymers, 
where failure occurs at applied strains of the order of 1%. 

Recently, an extended theory of dynamic fracture - the 
“Weakly Nonlinear Theory of Dynamic Fracture” - was 
developed [MS]. This theory still treats the large scales of 
a fracture problem as linear elastic, but explicitly takes 
into account the leading order nonlinear elastic correc¬ 
tions near the tip of a crack. In the latter region, strains 
are necessarily large and nonlinearities of the elastic fields 
are important for understanding the physical state that 
actually drives material failure. 

The weakly nonlinear theory introduced a new intrinsic 


length scale (i.e. independent of the external geometry 
and crack’s dimensions) that emerges from the compe¬ 
tition between linear and weakly nonlinear deformation. 
The new length scale, which is missing in LEFM and 
characterizes the scale near the crack tip where LEFM 
breaks down, has been argued to play an important role 
in determining the crack’s path [9] and has been shown 
to play a decisive role in crack instabilities [TOHI^ . This 
theory has recently received significant experimental sup¬ 
port [nmii. 

Of late, classes of materials in which large deforma¬ 
tions precede failure have become increasingly important 
and have been extensively studied in various contexts. 
These highly-deformable materials include synthetic elas¬ 
tomers, biopolymers, gels and a broad range of soft bi¬ 
ological tissues (e.g. arterial walls, veins, skin, tendons 
etc.). Reviewing even a small portion of the work done 
on the fracture of soft, highly-deformable, materials — 
which encompass various scientific disciplines and com¬ 
munities — goes well beyond the scope of this paper. We 
can only mention here a few examples. 

The fracture resistance of some synthetic materials, 
such as double-network hydrogels made of ionically and 
covalently crosslinked networks [THl420j . can be made 
to be enormous and new applications are expected to 
abound. Basic phenomena such as delayed fracture of 
soft solids na, and surface tension-induced and capillary 
fracture of gels [22U^ have been studied. Finally, the 
work of Baumberger and coworkers on the quasi-static 
fracture of soft materials has both elucidated the fracture 
mechanisms of biopolymer gels and established 

the existence of an intrinsic length scale associated with 
elastic nonlinearities HDHH]. 

It is important to note that many of the materials of 
interest here, although highly-deformable and tough, are 
still brittle. The materials that we are considering here 
are materials where high-deformability is dominated by 
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nonlinear elastic behavior, with little ductility. This type 
of behavior is typical in many of the new tough elas¬ 
tomers that are being developed. Since the bulk dissipa¬ 
tion can be ignored, these materials are indeed brittle - 
as the only dissipation takes place in the near vicinity of 
the crack tip. 

When failure occurs under the application of large 
background deformation (pre-stress/strain), the defor¬ 
mation is nonlinear everywhere in the material, invali¬ 
dating the assumption of remote linear elastic deforma¬ 
tion. To address this problem, some works considered 
the fully nonlinear field equations for certain classes of 
highly-deformable materials and non-perturbatively de¬ 
rived the leading asymptotic fields in the inner most crack 
tip region. The vast majority of these works focussed on 
static cracks though not all Such an 

asymptotic approach would be useful if the asymptotic 
fields are generic and universally linked to the remote 
loading, something which has not yet been established. 
Marder, following a different approach combining numer¬ 
ical and analytic techniques, developed a rather compre¬ 
hensive dynamic fracture theory of rubber m- Despite 
these important efforts, we are still far from having a 
well-established general theory of the dynamic fracture 
of highly-deformable, strongly nonlinear, materials. 

In this paper, we develop a theoretical framework in 
which the background, possibly finite (nonlinear), defor¬ 
mation induced by the external loading in the absence of 
a crack (pre-stress/strain) is treated non-perturbatively. 
Then, the effects of the crack on its near-tip region are 
treated perturbatively to second order with respect to the 
background deformation. This theory will be shown to 
offer quantitatively good approximations to experimental 
data at moderately large background deformation. Con¬ 
sequently, we occasionally refer to it as the moderately 
large deformation theory. 

The development of a perturbative fracture theory in 
the presence of non-infinitesimal background deformation 
is conceptually non-trivial. It raises the following ques¬ 
tion: if the background deformation is finite (nonlinear) 
and a crack significantly amplifies the background (re¬ 
motely applied) deformation near its tip, is it justified 
to treat the latter as a relatively small correction to the 
former? 

To quantitatively address this issue, we first develop 
the moderately large deformation theory and derive its 
near crack tip solutions (both analytically and semi- 
analytically). We then systematically and quantitatively 
compare it to the weakly nonlinear theory. Both theories 
are also compared to extensive direct measurements of 
the near-tip deformation fields of dynamic cracks propa¬ 
gating in a brittle elastomer gel, where the background 
deformation is increased in a controlled manner. 

Our results show that while the two theories agree 
with each other at relatively small background defor¬ 
mation, the moderately large deformation theory offers 
significantly better approximations to the experimental 
data at moderately large levels of external loadings (pre¬ 


stress/strain) and high propagation velocities. 

Theoretical frameworks that invoke perturbations of a 
pre-stressed state - which are sometimes referred to as 
“Mechanics of incremental deformation” - are not new 
in themselves (SSJ [31] . A classical example is the theory 
of small amplitude waves in nonlinearly pre-stressed ma¬ 
terials, which has recently attracted renewed attention 
[40) . In this case, a perturbative approach relative to the 
pre-stressed/strained state of the material is fully justi¬ 
fied as the wave amplitude can remain small relative to 
the large background deformation. As explained above, 
this is not automatically the case in fracture dynamics 
where the background deformation is significantly am¬ 
plified near the crack tip and hence in principle may not 
be treated as a small perturbation. 

A number of authors have previously discussed such 
approaches to fracture [39l [41], where the effect of the 
crack was treated perturbatively to linear order with re¬ 
spect to the background deformation. As far as we can 
tell, however, these authors did not address at all the 
range of validity of the approach. In particular, their per¬ 
turbative approach was confined to linear order, which - 
as we show below - is insufficient since higher order ef¬ 
fects (in particular, weakly nonlinear effects with respect 
to the background deformation) play an important role. 

Furthermore, to the best of our knowledge, quantita¬ 
tive comparisons of theoretical predictions to detailed ex¬ 
perimental data, as we do here, have not previously been 
performed. We believe that the combined theoretical- 
experimental results presented in this paper offer a useful 
framework to quantitatively address the fracture proper¬ 
ties and dynamics of highly-deformable materials. 

II. THEORETICAL FRAMEWORK 

To lay down the theoretical grounds for the approach 
we propose, consider a dynamic crack propagating in a 
2D nonlinear elastic solid described by an energy func¬ 
tional U{F). The deformation gradient tensor F is de¬ 
fined as F{x,t) = Wx‘p{x,t), where the motion (p{x,t) is 
a continuous, differentiable and invertible mapping be¬ 
tween a reference (undeformed) configuration described 
by X and a deformed configuration described by x\ such 
that x' — ip(x,t)=x+u(x,t). u{x,t) is the displacement 
vector field. Linear momentum balance can be expressed 
in the reference configuration as [12] 

Vx-s = pip , (1) 

where s is the first Piola-Kirchhoff stress tensor, s = 
dpU{F)^ and p is the time-independent reference mass 
density. 

The crack is assumed to follow a straight trajectory 
and to propagate steadily at a velocity v along the pos¬ 
itive x-axis. Under symmetric tensile loading along the 
y-axis, that is under mode-I fracture conditions, the crack 
faces are being separated and hence are traction-free. 
These traction-free boundary conditions on the crack 
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faces can be expressed in the undeformed configuration 
as [42] 

Sxyir, 0 = ±7r) = Syy(r, 9 = ±7r) = 0 , (2) 

where (r, 0) is a polar coordinates system co-moving with 
the crack tip (r = 0 is the tip location and 0 = 0 is the 
propagation direction). 

The solution of Eq. Q , with the boundary conditions 
of Eq. and a constitutive relation s = dpU{F) for 
a general nonlinear energy functional U{F), is analyt¬ 
ically intractable. To make progress, some approxima¬ 
tions are invoked, most notably in situations in which 
fracture occurs under small background deformation. In 
such cases, the displacement gradient tensor H = 
is treated as small everywhere except for a small zone 
around the crack tip, and a perturbative approach is de¬ 
veloped. LEEM [IHSj and the weakly nonlinear theory of 
fracture IlHS] fall under this category. Highly-deformable 
materials, on the other hand, fail at large, nonlinear, 
background deformation and in principle the status of 
a perturbative approach is not clear. 

A. General formulation 

To develop our approach, we write the motion tf(x, y, t) 
as 

ipx{x,y,t) = XxX +Ux{x,y,t) , 

(py{x,y,t) = Xyy+ Uy{x,y,t) . ( 3 ) 

Here we decompose the total motion into a contribu¬ 
tion emerging from the external loading Xx,y in the ab¬ 
sence of a crack (the pre-stress/strain), which are the 
stretches in the x and y directions respectively, and into 
the effect of the crack quantified by lA{x,y,t). Eor sim¬ 
plicity, and to later allow direct comparison with ex¬ 
periments, we assume hereafter that Xx,y are constants. 
When Xx = Xy = 1, U{x,t) becomes the ordinary dis¬ 
placement field u{x,t). In this case, LEEM corresponds 
to the linear approximation in the displacement gradient 
H and the weakly nonlinear theory to the second order 
approximation in H (leading order nonlinearity) 

The decomposition in Eq. ([^ is sketched in Fig. 

We now consider situations where Xx,y deviate to a 
finite degree from unity. When the motion of Eq. ([^ 
is substituted in s = dFU(F), which is then substituted in 
Eq. 0, a nonlinear set of equations in U is obtained. In 
this paper we solve this problem perturbatively to second 
order in i.e. employ the expansion 

U (4) 

in the near tip vicinity of a dynamic crack, which is un¬ 
derstood as a displacement gradients expansion. 

The expansion in Eq. @ corresponds to an expansion 
of the stress of the form s ~ 8^°^ -|- -|- where 

simply corresponds to the pre-stress. Using the latter 



FIG. 1. A sketch that illustrates the decomposition in Eq. §. 
(a) An undeformed configuration is uniformly stretched in two 
directions, described by two principal stretches Xx,y>l, cor¬ 
responding to the first terms on the right-hand-side of Eq. §. 
The effect of the crack relative to this pre-stretched state is 
described by the field U{x, y, t). (b) The actual motion is de¬ 
scribed by (fi{x,y,t) which takes into account both the pres¬ 
ence of the crack and the external stretch. On the left, the 
undeformed configuration with an unopened crack is shown. 
On the right, the stretched configuration with an opened crack 
propagating at a velocity v is shown. The coordinate system 
showing the propagation direction {x), the tensile loading di¬ 
rection (y) and the thickness direction (z) is added for clarity. 

and transforming into a frame of reference moving with 
the crack tip, momentum balance in Eq. 0 translates 
into two linear problems which take the form 

, (5) 

^ py^dxxU^^^ , ( 6 ) 

and Eq. 0 gives rise to the corresponding boundary 
conditions on and is a vectorial linear 

differential operator, which depends on the nonlinear en¬ 
ergy functional U{F) and Xx,y, and has the form of an 
effective anisotropic linear elastic operator. Anisotropy 
here emerges due to nonlinearity in U{F) and the pos¬ 
sibly symmetry-breaking external stretches X^^y Since 
C[-] is a linear elastic operator, we expect the asymptotic 
near tip solution of Eq. 0 to give rise to the standard 
singularity ^ l/v^ at small r. This square root 

singularity implies a parabolic crack tip opening profile 

m- 

What do we expect to occur to second order? We 
are guided by the intuition gained by solving the weakly 
nonlinear problem for infinitesimal deformations IHH]. 
Hence, the second order problem, i.e. Eq. 0 , fea¬ 
tures the same linear operator C[-\ as in the first order 
problem, but also an effective body force corre¬ 

sponding to quadratic contributions emerging from the 
first order solution. In particular, we have ^ 
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~ 1/r^. The boundary condition which 
satisfies, emerging from a consistent expansion of Eq. 

features an effective surface force proportional to 
l/r. As in the weakly nonlinear theory derived for small 
background strains, this structure is expected to lead to 
~ l/r, i.e. to a stronger singularity than the lin¬ 
ear problem. is expected to modify the crack tip 

shape and to introduce a new lengthscale into the prob¬ 
lem [IHS], the length at which \- 

The elastic fields transport a finite amount of energy 
into the tip region, as quantified by the J-integral [5] 


J = y {u{F) + \p[dt'^if'^vnx +SijTijdt^Pi dC , (7) 


where C is a contour encircling the tip and n is an out¬ 
ward unit vector on C. This integral is path-independent 
for steady-state crack propagation and for any contour C 
within a non-dissipative region described constitutively 
by the elastic energy functional U{F). G{v) = J/u, the 
energy release rate whose dimensions are energy per unit 
crack area, is dissipated near the tip. This dissipation is 
quantified by the fracture energy r(u) - a measure of the 
material’s resistance to crack propagation - which is a 
fundamental material function assumed to depend only 
on the crack propagation velocity. Energy balance im¬ 
plies that G{v) =r(z;), which enables us to use G{v) to 
calculate r(r;) (see below). 

The approximate solution for the motion <^(x, y) de¬ 
pends in a nontrivial way on the nonlinear energy func¬ 
tional U{F) and on the background stretches \x,y In 
the next subsections we will demonstrate how to actu¬ 
ally derive the solution and explore some of its physical 
properties. 


while the traction-free boundary conditions of Eq. ([^ 
read 


Sxy{r,6 = ±TT) = nidyifx + dxipy)\9=±n = 0 , 

^yy (^7^ — irr) — yi^dyipy dxPx ) |^—±77 — 0 ■ (IH) 

We focus now on pre-stressing/straining of the form 
Aa; = Ay = A, i.e. on a biaxial stretch, which is sim¬ 
pler than uniaxial stretching due to its symmetry. This 
loading corresponds to a pre-stress of the form = 

(A — Substituting Eqs. in Eqs. we 

obtain to linear order (again in the co-moving frame of 
reference) 


-b 3A-VV(V • = pv'^dxxU^^^ , (11) 


which has the structure Eq. ([^. Using Eqs. (10) we 
obtain the following boundary conditions at 0 = ±7r 

pdyUi^'^ +\-^yidxU‘i^^ = 0 , 

(1 -b Z\-^)pLdyU‘^^'> + 2X-^pdxUi^'‘ = 0 . (12) 


Obviously, LEEM is recovered in the limit A —>■ 1 (as the 
material is incompressible, the resulting LEEM plane- 
stress problem involves a Poisson’s ratio of ^). Equation 


(11) has the structure of an isotropic linear elastic prob¬ 


lem (Lame equation) with a A-independent shear modu¬ 
lus and a first Lame coefficient of the form (3A“® —l)/r. 
Consequently, we can employ rather standard complex 
functions techniques mm to obtain the following asymp¬ 
totic (near tip) analytic solution 


B. Analytic example of first order asymptotic fields 

To see how all of this works, we discuss an explicit ex¬ 
ample that can be worked out analytically in a rather 
straightforward manner. We consider an incompress¬ 
ible neo-Hookean material under plane-stress conditions, 
whose nonlinear energy functional takes the form |43j 

U{F) = ^ - 3) . (8) 

Here p, is the shear modulus, F is the 2D deforma¬ 
tion gradient and Az{x,y) = [det(F)]“^ = [dxPx dyPy — 
dxPy dyipx\~^ is the out-of-plane stretch. We choose this 
energy functional because it is relevant for many highly- 
deformable materials and it will allow us later to compare 
our predictions to direct experimental measurements. 

The stress tensor s corresponding to U{F) in Eq. ([^ 
reads s^- = p {^djPi — A^ Cik eji dipk), where e^- is the 2D 
alternator (i.e. exx = ^yy = 0, txy = —^yx = !)• The 
momentum balance of Eq. ([^ takes the form 

“f M ^X^y ^xA^ = PPx 1 

pSJ (py -b p dyPx ^y^z ^3:^3:] ~ P^y j (9) 
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U^^\r,e-,v,X) 


[(A ®+a^)v^cos(|^) - (1 +A ^)adas^sCos{^)] 

fiy/27TD[v^ A) 

+ A-«)V^sin (^)] . 


(13) 


The quantities ls,d and Og d are analogous to their 
standard LEFM counterparts, rather with a pre-stretch 
dependent dilatational wave-speed Cd(A) = \/l + 3A“®Cs, 
where the shear wave-speed Cg = yjiijp is unaffected 
by A. In particular, a.g^d = \J\-{vjCg^dY, ^g^d = 
— (usin0/cs_d)^ and tan(0s_d) = as,dtan(0). The pre¬ 
stretch dependent analog of the Rayleigh function takes 
the form I?(u, A) = 2(1 -|- X~^)adOis ~ (l + Q^s) + <^ 5 ) 
and Kj is the mode-I stress-intensity-factor ^ [S]. The 
solution in Eq. (13) features the expected singularity, 
'^l/y/r, where the standard LEFM solution [5] is 
recovered in the limit A —>■ 1. 


The analytic solution in Eq. (13) has several phys¬ 


ical implications. First, it can be used to calculate the 
J-integral of Eq. ([^, employing the linear elastic approx¬ 
imation of U{F), yielding 

(14) 

The stress-intensity-factor Kj cannot be obtained from 
the asymptotic solution, rather from the global boundary 
value problem. It can be calculated analytically only in 
relatively simple cases and in general it is obtained nu¬ 
merically or measured experimentally. Once it is avail¬ 
able, Eq. ( |I4| ) allows one to calculate the fracture en¬ 
ergy r(u), a basic material property, through the relation 
T{v) = G{v). 

Alternatively, if r(w) is known (either from a proper 
dissipation theory, which is very rare, or through inde¬ 
pendent measurements) one can calculate Kj{v, A) using 
energy balance and Eq. (14). This clearly demonstrates 


that Kj depends on A; not being aware of this pre-stretch 
dependence can induce mistakes. 

Finally, Eqs. (§, (|3 and ( [T^ can be used to calculate 
the shape of the tip (often called crack tip opening dis¬ 
placement/profile) as ipx{r,Tr) = —k{v, X)(py{r,Tr), where 
the tip curvature reads 


i{v,X) = 


fiy/^D{v, A) 

2ad(u, A)(a2 - 1) 


iff 


(15) 


Since the tip curvature k is, in principle, a directly mea¬ 
surable quantity, the last result can be used to extract 
the stress-intensity-factor Kj. Again, we see that not 
being aware of the A-dependence (e.g. using instead the 
A = 1 result) will lead to mistakes. More generally, anal¬ 
ysis of the solution in Eqs. (13)-(15) reveals that the 


pre-stretch A has a marked effect on various important 
physical quantities, and that this effect increases signifi¬ 
cantly with increasing propagation velocity v. 


Up to now we considered the asymptotic solution to 
linear order in lA, which provided us with some insight 
into what kind of effects can be associated with the 
pre-stretch (related points were made in [39l ST]). We 
know, however, that even for small background deforma¬ 
tion second order nonlinearities are essential SHE). Con¬ 
sequently, in the next subsection we consider the solu¬ 
tion for both and under uniaxial background 
stretch. 


C. Semi-analytic second order asymptotic solution 

Here we focus on pre-stressing/straining of the form 
Xx= A“^/^ and Ay = A, i.e. on uniaxial stretch, which is a 
more commonly used experimental loading configuration 
(note that the out-of-plane pre-stretch is A 2 = A“^/^, en¬ 
suring incompressibility, XxXyXz = !)• It corresponds to a 
pre-stress Syy^ = (A—A“^)/r, where the other components 
of 5 *^°) vanish. Following the same procedure as above, 
Eq. (§ takes the form 


idxxHi^^+SX-^^^dxylA^^^ +dyylAi^^ = "^dxxUi^^ , ( 16 ) 

with the following boundary conditions at 0 = ±7r 

+ X~^/'^^dxUy^^ = 0 , 

(1 -h iX-^)^l^yU'^^'^ + 2X-'^l'^[idxU^J^ = 0 . (17) 


Equations @-([17])’ which manifestly exhibit elastic 
anisotropy, can in principle be solved analytically in the 
asymptotic regime of small r. The solution, however, 
is rather lengthy and we present here instead a semi- 
analytic procedure to obtain it. To leading order in small 
r, we expect ^ y/r and the angular dependence to 
be expressed as a half-integer Fourier series. The sub¬ 
leading term in small r (i.e. in the expansion in space), 
which has not been discussed up to now and which will 
be included below to enable direct comparison with the 
experiments to follow, makes a contribution oc r cos 9 to 
lAx^^ and oc rsin 0 to lAy^^ (the boundary conditions de¬ 
termine the ratio between the amplitudes). 
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Therefore, we have 


uW{r,9)=Uo + KiJ—^< 

V ^ n=l 
T 


, COS 


(2n- 1)9 


H-(1 + 3A "^Ircos^ 

12/i 


U^^^\r,9)=Ki. 


sin 

■(2n- 1)9' 


2 


- —X-^/'^rsm9 

6fi 


(18) 


Here Uq is a constant does not include such a con¬ 

stant due to the mode-I symmetry), Kj is a dimensionless 
stress-intensity-factor and {a„,6„} are dimensionless co¬ 
efficients. The coefhcients of the sub-leading term have 
been chosen so as to satisfy the boundary conditions of 
Eq. (17). Note that in the limit A —?> 1 this term cor¬ 


responds to the so-called T-stress and consequently we 
have chosen the yet undetermined amplitude to agree 
with the standard result in this limit, where T is a quan¬ 
tity of stress dimensions m- 

The coefficients {a„, are determined by the 

set of linear algebraic equations obtained upon substitu¬ 
tion of Eqs. (18) in Eqs. (|I^-([T7|), where N is chosen to 
be sufficiently large to ensure convergence. Kj is deter¬ 
mined by evaluating the J-integral in Eq. ([^ (with the 
linear elastic approximation of U{F)) and equating G{v) 
to the fracture energy r(u). T will be extracted from 
experimental data. 

As stated above, previous work has conclusively 
demonstrated that second order nonlinearities are im¬ 
portant [MS]. Hence we wish to calculate To that 
aim, we follow the procedure described above to calcu¬ 
late from which Eq. can be obtained. A very 
detailed, step-by-step, explanation of the mathematical 
procedure can be found in section 4.2 of [8] and in [9]. 
Based on the solution obtained in the framework of the 
weakly nonlinear theory of fracture SHE], we expect the 
solution for to take the form 

(a ^^0 log r -t ^ c„ cos [n6»]^ , 

U^^\r,9) = l{^o9+ ^driSi\\[n9^ . (19) 

Note that, in principle, the argument of log r should have 
been made non-dimensional, but this would simply rede¬ 
fine Uq in Eq. (18) and hence is not essential. 

The solution in Eq. (19) has the property that 
^ 1/r. As was shown previously, this singular¬ 
ity is special in the sense that it can produce a spurious 
force in the crack parallel direction lilllll]. To eliminate 
it, we supplement the equations of motion and boundary 
conditions with the additional constraint 


fx = 


s^^}{r,9) cos 9 + 9) sin 9 rd9 = 0, (20) 


where fx is the net force per unit sample thickness acting 
in the x direction on a line of radius r encircling the crack 

tip [HI [7]. 

By satisfying all of these equations, we can calculate 
the coefficients {cn,dn}n=o-N, where N is chosen to be 
sufficiently large to ensure convergence. Note that the 
equation for cf. Eq. (|^, as well as the boundary 

conditions, require the knowledge of Once done, 

the solution in Eq. in the asymptotic region of small 
r, is at hand. In the limit A —> 1, the weakly nonlinear 
theory of fracture is recovered. The next step will be to 
quantitatively test the predictions of the theory devel¬ 
oped above against direct experimental measurements at 
various values of A and v. 


III. COMPARISON TO EXPERIMENTS 

Our goal here is to experimentally test the new theo¬ 
retical framework. In particular, we would like to com¬ 
pare the predictions of the weakly nonlinear theory (ob¬ 
tained in the limit A —>■ 1 in the theory above) to the 
predictions of the moderately large deformation theory 
(which is formulated relative to a finite A), against ex¬ 
perimental data. To that aim, we conducted experiments 
with polyacrylamide gels, which are transparent, homo¬ 
geneous, brittle and incompressible elastomers. The gel 
composition used here is 14% aerylamide/bis-acrylamide 
with a 2.7% cross-linker concentration, providing a shear 
modulus fi = 32.3 ± 1.6KPa and a shear wave-speed 
Cs = 5.6 ± O.I5m/sec. /r is measured prior to each exper¬ 
iment to mitigate any small variations of the gel prop¬ 
erties. The energy functional in Eq. (|^ quantitatively 
describes the gel [55] . 

The typical dimensions of our samples are {x x y x z) 
120 X 120 X 0.3mm in the crack propagation, tensile load¬ 
ing and thickness directions, respectively. The thick¬ 
ness was chosen to statistically suppress micro-branching 
[45] : while micro-branching can occur for all velocities 
0.3cs < V < 0.9cs, the probability of exciting them de¬ 
creases for both thin samples and increased crack accel¬ 
erations. All of the results presented are for single-crack 
states in which micro-branching is not observed. 

The experiments are performed under uniaxial tensile 
loading in the y-direction in accordance with the uniax¬ 
ial theory presented in the last section. The sample was 
held at a constant stretch in the range A= 1.058 — 1.129 
prior to crack initiation at the mid-edge of the sample’s 
vertical boundary, as described in Fig. Measurements 
of the crack and its surrounding displacement fields were 
made with a fast camera (IDT-Y7) focused on an area 
of 17.4 X 9.8mm with a spatial resolution of 1920 x 1080 
pixels and frame rates of 7400 — 8100 frames/sec. The 
crack velocity for each measurement was set essentially 
by varying both the imposed strain and position of the 
measurement area used. All of the cracks are acceler¬ 
ating in accordance with the initial strain imposed by 
displacing the vertical boundaries of the sample. We 
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achieved desired velocity ranges in each experiment by 
varying the location of the measurement area according 
to the initially imposed strain. Due to the small size of 
the measurement area, the crack velocities were approx¬ 
imately constant throughout the measurement region. 

As in [141 Efi] j the gels are cast in a mold upon which 
a rectangular grid is printed on one of its xy surfaces. 
The grid was formed by lithographic printing on a spin- 
coated epoxy layer. This process created a perfect square 
mesh of depth 2/rm (in the z-direction) and lattice spac¬ 
ing 60/im (in the xy plane). Upon casting, this grid 
is imprinted on one face of the gel sheets as shown in 
Fig. Shadowgraphy, using strobed lighting (2/isec 
duration), is used to image both the deformed grid and 
crack opening profile as a crack’s tip propagates across 
the field of view (Fig. The location of the center 

of each grid point in the deformed grid is determined to 
within lO^m resolution. The displacement fields were ac¬ 
quired by comparing the position of the grid points under 
deformation to their position in a deformation-free sys¬ 
tem. Fig. demonstrates a typical measurement of the 
displacement-gradient field (here the Syy = dyUy compo¬ 
nent is shown). 



FIG. 2. (a) Experiments are performed with transparent thin 
sheets of a brittle elastomer gel with an imprinted grid on one 
xy face {y is the tensile loading direction and x is the crack 
propagation direction). Collimated light passed through the 
sample enables shadowgraph visualization of the grid while 
pictures are taken with a fast camera, (b) Typical photograph 
of a crack propagating at u = 0.53cs through the grid, (c) The 
measured displacement-gradient field component eyy{x,y) = 
dyUy{x,y) extracted from panel (b). 

In Fig. [^we compare the experimental measurements 
with both the weakly nonlinear theory and the moder¬ 
ately large deformation theory for a moderate stretch 
A = 1.058 and crack propagation velocity v = 0.29cs. The 


basic field of interest is the displacement field u{x,y), 
which is measured directly and is theoretically obtained 
from Eq. ffl through Ux{x,y) = {X~^^'^ — 1) x +Uy;{x,y) 
and Uy(x, ^={X— l)y +Uy{x, y). Each theory contains 
three parameters {Uq,T,T} that are not determined by 
the asymptotic analysis. U^s and T correspond, respec¬ 
tively, to a small constant shift of the crack tip location 
and to the T-stress. As our measurements are not solely 
within the asymptotic (singular) region, both of these 
quantities are needed for a good quantitative compari¬ 
son. 

For both theories, we determined the values of 
{Wo,r, F} by the following procedure. We considered 
both Ux{r = x,6 — 9) and 0 = Tr) [ipy (r, 9 — ±7r)] (in 

the latter, r parameterizes the function), i.e. the crack 
parallel displacement ahead of the tip and the crack tip 
opening profile, respectively. The first function is quite 
sensitive to T, whereas the latter is very sensitive to F. 
We iteratively performed a 3-parameter fit over these two 
functions until the best fit with the same T and F is ob¬ 
tained. Uq is not constrained to be the same, but turns 
out to be so. The resulting fits for Ux{r = x,9 — 0) 
and (fx{r,9 = tt) [ipy{r,6 — ±TT)] are shown in Figs, 
b, respectively. The convergence in the T — T paramet¬ 
ric plane, for both the weakly nonlinear and moderately 
large deformation theories, is demonstrated in the inset 
of Fig. [^. 

Figures |^-b suggest that at this level of pre¬ 
stressing/straining (relatively low in the context of the 
results to follow) and crack propagation velocity the two 
theories appear to be almost indistinguishable, at least 
as far as the quantities shown are considered. The pa¬ 
rameters, however, are not the same. In particular, as 
the inset of Fig. clearly demonstrates, in this case F is 
quite similar while T is not. The latter difference is ex¬ 
pected since the background stretches in the moderately 
large deformation theory, cf. Eq. (|^, directly affect the 
T-stress term and actually make T positive (while it is 
negative in the weakly nonlinear theory). We emphasize, 
though, that when considered relative to the undeformed 
configuration, the T-stress in the moderately large defor¬ 
mation theory is also negative (as is common for uniaxial 
loading). 

Once {77o,T, F} are determined, there are no longer 
any free parameters and the two theoretical frameworks 
can be independently tested against other experimentally 
measured quantities. In particular, we will use two types 
of tests: 

• Comparing the predictions of the two theories to the 
measured tensile strain ahead of the tip, £yy{x,y = 9) = 
dyUy{x, y)\y = 0- 

• Comparing the predictions of the two theories for the 
fracture energy to r(?;) = G(v) = independently ob¬ 
tained from the J-integral in Eq. Q (using the measured 
ip{x,y), as in [3n])- 

The results of these parameter-free comparisons are 
shown in Fig. (main panel) and Fig. (focus on 
?; = 0.29cs). Figure]^ shows that the predictions of the 
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FIG. 3. Measurements of a crack propagating at v = 0.29cs 
under uniaxial stretch of A = 1.058. (a) The measured dis¬ 
placement in the crack parallel direction ahead of the tip 
Ux{r = x,6 = 0) (blue circles). Fits to the weakly nonlinear 
theory (red solid line) and to the moderately large deforma¬ 
tion theory (black dashed line) are superimposed, (b) The 
measured crack tip opening profile (fix{r,9 = n) [ipy(r,6 = ±.'K)] 
(blue solid line). As in panel (a), the fits to the two theories 
are superimposed. The fitting parameters {Uo,T, F} were ob¬ 
tained for each theory, where the convergence of the iterative 
procedure in the F — T parametric plane is shown in the inset 
of panel (c), for the weakly nonlinear theory (left) and mod¬ 
erately large deformation theory (right). The intersections in 
the F — r plane, denoted by black squares, are the values 
chosen by our fitting procedure (see text). In addition, we 
obtained Wo — 25/rm for both theories. Note that the green 
and blue curves were obtained by the hts to the crack open¬ 
ing profile and Ux(x,y = G), respectively, (c) The measured 
tensile strain ahead of the tip eyy{x,y = 0) = dyUy(x,y)\y=o 
(blue circles). The predictions of the two theories, using the 
parameters obtained by the fits displayed in panels (a) and 
(b), are superimposed (lines as above). 


two theories are in good agreement with the measured 
data, where the weakly nonlinear theory is doing slightly 
better. Figure shows that for v = 0.29cs, the fracture 
energies F predicted by the two theories are in agree¬ 
ment with the one independently calculated through the 
J-integral. All in all, we conclude that for relatively low 
levels of pre-stressing/straining and crack propagation 
velocities, the two theories appear rather consistent with 
one another and quantitatively agree with the experi¬ 
ments. Furthermore, this analysis reconfirms the valid¬ 
ity of the weakly nonlinear theory, as reported previously 
mi]. The main question now is what happens as the 
pre-stress/strain and the crack propagation velocity are 
significantly increased. 

In Fig. 1^ we repeat the analysis presented in Fig. 

for a crack propagating at v = 0.74cs under signifi¬ 
cantly increased pre-stressing/straining corresponding to 
A = 1.096. Figures & -b indicate that both theories can 
be reasonably fitted to the measured Ux{r = x,9 = 0) and 
‘Px 0 = tt) [py (r, 9 = ±7r)] , where the moderately large 
theory is doing better with respect to the former. The 


FIG. 4. Measurements of a crack propagating at u = 0.74cs un¬ 
der uniaxial strain of A = 1.096, where everything is as in Fig. 
1^ except for the inset of panel (c). In the inset, we plot the 
difference between predictions of the moderately large defor¬ 
mation and weakly nonlinear theories for eyy{x, y = 0) (shown 
in the main panel), denoted by Asyy. The difference is a non¬ 
trivial spatially-varying function that increases significantly 
as the crack tip is approached. 

parameter-free comparison shown in Fig. |^, however, re¬ 
veals a striking difference between the two theories; the 
prediction of the moderately large deformation theory for 
£yy{x, y = 0) is significantly better than the prediction of 
the weakly nonlinear theory and is in good quantitatively 
agreement with the measurements. Furthermore, Fig. 
shows that the moderately large deformation theory pre¬ 
dicts a fracture energy F much closer to the independent 
J-integral estimate than the weakly nonlinear theory (fo¬ 
cus on u = 0.74cs). 

It is crucial to understand that the pre-stretching A af¬ 
fects all physically relevant quantities in the problem, as 
in the analytic example of Sec. |II B[ in a nontrivial way. 
For example, in the inset of Fig. we show the differ¬ 
ence between the predictions of the two theories shown 
in the main panel, which is a nontrivial spatially-varying 
function that increases significantly as the crack tip is 
approached. The quantitative analysis presented in Figs. 

and has been repeated for many cracks with prop¬ 
agation velocities in the range v = 0.26cs — 0.75cs and 
pre-stretching levels in the range A = 1.058— 1.129 (since 
the cracks are mildly accelerating, each pre-stretch pro¬ 
duces a range of crack velocities). As in Fig. |^, at 
stretches of about A = 1.1 or higher and large propaga¬ 
tion velocities, the moderately large deformation theory 
predicts £yy{x,y = 0) significantly better than the weakly 
nonlinear theory. 

Before we discuss the predictions for fracture energy 
F(u), we briefly highlight the importance of second or¬ 
der nonlinearities in the expansion relative to the pre¬ 
stretched configuration. This has been previously estab¬ 
lished in relation to the weakly nonlinear theory, cf. Fig. 
1 in [3]. To show that this remains valid in the case of 
the moderately large deformation theory, we plot in Fig.[^ 
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FIG. 5. The same as Fig. except that the weakly nonlin¬ 
ear analysis is replaced with a moderately large deformation 
analysis truncated to first order (green solid line). For the 
latter, we follow the same fitting procedure as before, but 
set = 0. The best htting parameters are Uo = 100/rm, 
r/3/i = 0.01 and F = 19.8 J/m^. The necessity of the sec¬ 
ond order nonlinearities is evident from the panels (b) and 
(c), where signihcant discrepancies between the first order 
truncated moderately large deformation theory and the ex¬ 
perimental data are observed. 
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FIG. 6. The measured fracture energy F(ti) (squares) with 
propagation velocities in the range v = 0.26cs — 0.75cs and 
pre-stretching levels in the range A = 1.058 — 1.129 (different 
colors correspond to different A, cf. legend). r(?;) was es¬ 
timated from the J-integral of Eq. using the measured 
(p{x,y). The predictions of the weakly nonlinear theory (tri¬ 
angles) and of the moderately large theory (circles) are pre¬ 
sented for comparison. 
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everything as in Fig|§ except that the weakly nonlinear 
analysis is replaced with a moderately large deformation 
analysis truncated to first order. For the latter, we fol¬ 
low the same htting procedure for {Uq.T.T} as before, 
but set = 0. The results clearly demonstrate that 
second order nonlinearities are indeed essential. In par¬ 
ticular, the footprints of the missing logarithmic term 
(which appears in Ux , cf. Eq. (19)) and the stronger 
singularity ~ 1/r) are evident in Figs. [lo-c. The 

rather dramatic failure of the linear order moderately 
large deformation theory to predict duUy{x,y = 0) is re¬ 
lated to the properties of the LEFM asymptotic helds, 
which predict a negative duUy{x,y = 0) for sufficiently 
high velocities [3]. 

Figure]^ shows r(u) for the full range of pre-stretches 
and propagation velocities considered in this study, in¬ 
dicating that the moderately large deformation theory 
predicts fracture energy values closer to the independent 
J-integral estimate than the weakly nonlinear theory (ac¬ 
curately predicting the J-integral values up to ^ 10%). 
We would like to stress that the performed fits are ro¬ 
bust. In fact, we checked that using eyy{x,y = 0) and 
— tt) [ipy{r,0 = ±TT)] to determine {lA[j,T,T} and 
then testing the predictions for Ux{r = x,0 — 0) and F 
(independently obtained through the J-integral) yielded 
similar results to those reported above. 

Taken together, we believe that the quantitative analy¬ 
sis presented here clearly shows that the moderately large 
deformation framework offers significantly better approx¬ 
imations to the direct measurements than the weakly 
nonlinear theory at moderately large levels of external 
loadings and high propagation velocities. 


IV. CONCLUDING REMARKS 

In this paper we developed and experimentally tested 
a dynamic fracture theory of highly-deformable mate¬ 
rials which fail under the application of large external 
strains. The theory is based on a second order expan¬ 
sion in the displacement-gradients with respect to a non- 
linearly stretched reference state. While the theory is 
mathematically well-defined, its physical range of appli¬ 
cability - and the mere existence of such a range - are 
not a-priori guaranteed. The reason for this is that as 
the background strains are already rather large, ampli¬ 
fication of deformation near the crack tip may render a 
perturbative approach inappropriate. 

Our direct experiments showed firstly that the back¬ 
ground strain has a significant effect on the deformation 
fields surrounding crack tips. Even at moderate strains, 
the standard theory fails to provide a good description of 
the near tip fields. Secondly, and quite surprisingly, the 
experiments showed that the new theory provides a good 
description of these fields. This central finding emerged 
because the amplihcation of deformation in the region of 
interest was not enormous and due to the inclusion of 
second order terms in the theory. These results imply 
that the theory may have a robust range of applicabil¬ 
ity, at least at moderately large background strains and 
high propagation velocities, but possibly also at larger 
background deformation. 

This convincing experimental support indicates that 
the proposed theoretical development offers a framework 
to understand the dynamic fracture of soft materials that 
fail under large pre-stressing/straining (of the order of 















10 


~ 10—20% or larger), going significantly beyond the stan¬ 
dard fracture theory of ordinary materials that fail under 
strains of ~ 1%. The theory shows that in order to quan¬ 
titatively understand the helds that drive material fail¬ 
ure near crack tips, the deformation-induced anisotropy 
and fundamental material properties such as the fracture 
energy, the pre-stressing/straining needs to be properly 
taken into account in such materials. The theory may 
find applications in a range of problems dealing with the 
failure of soft materials, from food processing to tissue 
rupture. 

It is important to note that while some existing liter¬ 
ature exclusively focusses on the inner most asymptotic 
crack tip region in highly-deformable materials P51l35j . 
our theory takes into account the pre-strained/stressed 
large scales and links them to the near tip region as ap¬ 
proached from the outside. As such, the theory should 
be regarded as intermediate asymptotics. The experi¬ 
ments presented above, which quantitatively support the 
theory, are able to probe this intermediate asymptotic 
region. 

One insight emerging from this work is that the form 
of these intermediate asymptotic solutions is highly influ¬ 


enced by the magnitude of the background deformations. 
In the classic LEFM theory, the loading and background 
stresses are solely accounted for by their influence on the 
intensity of the universal singularity, the stress intensity 
factor. Here, we have shown that large strains inherent in 
the external loading actually influence the fields in a va¬ 
riety of subtle ways that can not solely be accommodated 
by a change in the value of the stress intensity factor. 

Finally, note that various highly-deformable materials 
fail under background deformation of the order of 100% 
or even larger, a regime that has not been probed by 
the experiments presented above. Therefore, it remains 
to be seen in future work how far one can push the 
theoretical framework developed here. 
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